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Abstract 

There has been much research on codes over Z 4 , sometimes called 
quaternary codes, for over a decade. Yet, no database is available for 
best known quaternary codes. This work introduces a new database 
for quaternary codes. It also presents a new search algorithm called 
genetic code search (GCS), as well as new quaternary codes obtained 
by existing and new search methods. 
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I Introduction 

One of the main problems of coding theory is to construct codes with best 
possible parameters. There are databases of best known codes over small 
finite fields. For many years the online table [9] has been the primary source 
of the records of the best known codes over small fields. Recently, it is 
announced that this table is discontinued due to the existence of |15j which 
often has more explicit information on constructions. The computer algebra 
system MAGMA [ 8 ] has such a database too. Moreover, a table of best 
known binary non-linear codes is available at |18j . 

For over a decade there has been intensive research on codes over Z 4 , 
integers modulo 4, sometimes called quaternary codes. The term “quater¬ 
nary code” has been used both for codes over Z 4 and for codes over F 4 , the 
finite field with 4 elements. In this paper, we shall use the term exclusively 
for Z 4 codes. Among other results, some good quaternary codes have been 
constructed [TO] , [7] , [23] , and [5]. Self-dual codes over Z 4 of length up to 9 
are classified in and this is extended to length 15 in [14] (16 for Type 

II codes in [20]). Rains has classified optimal self-dual codes over Z 4 in 
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|22j . A large number of self-orthogonal quasi-twisted (QT) Z 4 codes have 
also been constructed m- Despite all this research, no database of best 
known quaternary codes is available. The development of such a table has 
been started in [5]. We now have compiled a database of quaternary codes. 
It is available at http;//Z4Codes.info/ and it is being updated continually. 
Unlike the tables at m, we do not overwrite the existing entries when they 
are improved but rather keep both the old and the new results. This strat¬ 
egy is chosen primarily based on the fact that several different metrics on 
quaternary codes, i.e. Hamming, Euclidean and Lee distance, have been 
considered by the researchers. Moreover, for the sake of easier communica¬ 
tion, we have decided to provide researchers with administrative privileges 
that would allow them to add their new results to the table, as well as edit 
the existing entries. Accounts can be acquired by contacting the database 
editors via email. 

In addition to the creation of the database, we have also devised and 
implemented some search algorithms to find new quaternary codes. The 
paper also describes and reports the results of these searches. One of the 
search methods we consider in this paper is the further exploration of the 
class of quasi-cyclic codes, which has been the source of many of the new 
codes discovered in recent years. The other method is called the “progressive 
dimension growth” (PDG) which is introduced recently in [3] for helds. In 
our work, it is modified for the ring Z 4 and the Lee metric. Finally, extending 
some of the ideas behind PDG, we implemented a new algorithm called 
“genetic code search” (GCS) that has produced better results than PDG 
over Z 4 . The following sections give more information about the algorithms. 
We have used MAGMA for all computations. 

2 Basic Facts on Quaternary Codes 

A code C of length n over Z 4 is a subset of ZJ. C is a linear code over Z 4 if 
it is an additive subgroup of Z 4 , hence a submodule of Z 4 . In this paper we 
will consider only linear codes over Z 4 . An element of C is called a codeword 
and a generator matrix is a matrix whose rows generate C. The Hamming 
weight wh{x) of a vector x = {xi,X2, ■ ■ ■ ,Xn) in ZJ is |{i : Xi 7 ^ 0}|. The 
Lee weight wl{x) of a vector x is niin{|xj|, |4 — Xi|}. 

The Hamming and Lee distances dnix, y) and dL{x, y) between two vec¬ 
tors X and y are wh{x — y) and wl{x — y), respectively. The minimum 
Hamming and Lee weights, djj and di, of C are the smallest Hamming and 
Lee weights, respectively, amongst all non-zero codewords of C. 
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The Gray map (j) : 1^2 ^ is the coordinate-wise extension of the 
function from Z 4 to dehned by 0 —>■ (0,0), 1 —)■ (1,0), 2 —> (1,1), 3 ^ 
(0,1). The image, 4>{C), under the Gray map of a linear code C over Z 4 
of length n is a (in general non-linear) binary code of length 2n. The Gray 
map is an isometry from {7i2, wl) to {'Z 2 "',wh)- Therefore, the minimum 
Hamming weight of 4>{C) is equal to the minimum Lee weight of C. 

Two codes are said to be equivalent if one can be obtained from the 
other by permuting the coordinates and (if necessary) changing the signs of 
certain coordinates. Codes differing by only a permutation of coordinates 
are called permutation-equivalent. Any linear code C over Z 4 is permutation- 
equivalent to a code with generator matrix G of the form 


Iki Ai Bi -\- 2B2 
0 2/^2 2A2 


( 1 ) 


where Ai, A 2 , Bi, and B 2 are matrices with entries 0 or 1 and 1^ is the 
identity matrix of order k. Such a code has size 4 ^ 1 2^^. The code is a free 
module if and only if /c 2 = 0. If C has length n and minimum Lee weight 
di, then it is referred to as an ^drYcode. 


2.1 Cyclic Codes over Z 4 

A cyclic code over Z 4 is a Z 4 -linear code which is invariant under cyclic shifts 
where the cyclic shift of an m-tuple (xq, xi,..., Xm-i) over Z 4 is the m-tuple 
(xm_i, xo, • • •, Xm- 2 )- Similarly to the case of finite helds, cyclic codes over 
Z 4 of length n are ideals in the ring under the usual identification 

of vectors with polynomials. Although algebraically cyclic codes have the 
same structure over fields and over Z 4 (ideals in a factor ring), the fact that 
Z 4 [x] is not a unique factorization domain makes it more challenging to find 
all cyclic codes over Z 4 . For instance, computer algebra systems (such as 
Magma and Maple), cannot directly provide factorizations of x” — 1 for an 
arbitrary n. When n is odd, it is easier to obtain a factorization of x” — 1 
and hence to find all cyclic codes of length re. For an even re, the situation 
is much harder. In fact, the factorization is not unique in that case. In this 
paper we consider only cyclic codes of odd length over Z 4 . For the case of 
odd re some of the most important facts about ideals of the relevant ring and 
the factorization of x"' — 1 are summarized below, and they can be found in 
[21], [ 25 ] or [27]. For the case of even re, we refer the reader to [I],[2], and 
i. 

For an odd positive integer re, x” — 1 can be factored into a product of 
finitely many pairwise coprime basic irreducible polynomials over Z 4 . Also, 
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this factorization is unique up to ordering of the factors 1211127!. In fact, we 
have the following: if f 2 {x)\{x"' — 1 ) in 1 j 2 [x] then there is a unique, monic 
polynomial f{x) G Z^[x] such that /(x)|(x” — 1 ) in 7,i[x] and f{x) = f 2 {x), 
where f{x) denotes the reduction of f{x) modulo 2 [27] . The polynomial 
f{x) is called the Hensel lift of / 2 (x). There are well-known methods of 
finding this polynomial, such as Graeffe’s method [27]. Therefore, there is 
a one-to-one correspondence between irreducible factors of x”' — 1 over Z 2 
and irreducible factors of x” — 1 over Z4. 

Once the factorization of x"" — 1 over Z 4 is obtained, the ideals of R := 
(x"-i} determined. For an odd positive integer n, any ideal I of 

the ring R has a generator of the form a / = {f{x)h{x),2f{x)g{x)) where 
f{x)g{x)h{x) = x" — 1 [lUllT]. Moreover, |/| = sG). it follows 
that the number of cyclic codes of length n is 3'’, where r is the number of 
irreducible factors of x” — 1 | 21 j . 

Finally, it can be shown that any ideal of R, for an odd n, is a principle 
ideal, with a generator of the form p(x) = f{x)h{x) + 2f{x) (or equivalently 
p{x) = f{x)h{x) + 2f{x)g{x)) where f{x),g{x),h{x) are as above [2X1 [27] . 

Remark 1 When x"" — 1 has r irreducible factors over a field, the total 
number of cyclic codes is 2^. We have a larger number over Z4 due to the 
existence of non-free codes (over a field all codes are free). 

Remark 2 The generator polynomial p{x) of an ideal of R described 
above does not necessarily divide x"’ — 1. For example, let n = 3, / (x) = 1, 
and h{x) = x — 1, then p{x) = x -|- 1 and p{x) /f(x^ — 1). When h{x) = 1, 
p{x) = 3/(x) = —/(x) does divide x” — 1. It is shown in [7] that the cyclic 
code generated by p{x) is a free module if and only if p{x) divides x” — I. 

2.2 Quasi-Cyclic Codes over Z4 

Much research has focused on the class of quasi-cyclic (QC) and the related 
class of quasi-twisted (QT) codes, and many new codes over small finite 
helds have been discovered within these classes. Some of these results can 
be found in mmm . [T 6 ] and [23]. In addition to the case of a field, QC 
codes over rings, especially over Z4, have been studied as well. QC codes 
over Z4 are first studied in [7], where a number of “good” quaternary codes 
are obtained. A quaternary linear code C with parameters [n,4^^2^^,d] 
(where d is the Lee weight) is called good ii d > d', where d' is the minimum 
distance of a best known binary linear code of length 2 n, and dimension 
2ki + k 2 , i.e. if the Gray image of C has a larger minimum distance than the 
comparable binary linear code. Similarly, a quaternary code will be called 
decent if its Gray image has the same parameters as the best known binary 
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code (i.e., if d = d'). The reason for this type of comparison is that even 
though the Gray image of a quaternary code is most likely non-linear, we 
do not have any other means of testing how good the parameters of a Z 4 
code are due to the facts that 

a) the table [TS] is much smaller than the tables for binary linear codes (it 
only goes up to minimum distance 29), and 

b) there are no extensive tables of quaternary codes. 

We believe that the table presented in this paper will meet a need in this 
area. The researchers are welcome to report and enter their codes to this 
database. 

Next, we summarize some of the basic facts concerning the structures of 
QC codes. A more detailed treatment can be found in [ 6 ] for QC codes over 
fields, and in [7] for QC codes over Z 4 . A linear code over a ring is called 
/—QC if it is invariant under the cyclic shift by I positions. Algebraically, an 
/-QC code of length n = ml over a ring R can be viewed as an R[x]/ (x™ — 1) 
submodule of {R[x]/{x"^ — 1))^ Then an r-generator QC code is spanned 
by r elements of {R[x]/{x^ — 1))^ In this paper, as is the case in most of 
the literature, we restrict ourselves to 1 -generator QC codes. The following 
is a generalization of an important result about 1 -generator QC codes [ 6 ], 
[7] that has been used in many of the recent work |12].|13|. The ring R can 
be a finite field, or Z 4 . 

Theorem 2.1 Let C be a 1-generator l-QC code of length n = ml with a 
generator of the form: 

9{x) = {fi{x)g{x), f 2 {x)g{x),fi{x)g{x)) ( 2 ) 

where g{x)\{x'^ - 1), g{x), fi{x) € - 1), and {fi{x),h{x)) = 1, 

h{x) = for all 1 < i < 1. Then I ■ d < d{C), where d is the minimum 

distance of the cyclic code generated by g[x), and d{C) is the minimum 
distance of C. Moreover, the dimension of C is equal to the dimension of 
the cyclic code generated by g{x). 

In terms of generator matrices, the QC codes can be characterized as 
follows. 

Let 



90 

91 

92 ■ 

9m —1 


9m —1 

90 

91 ■ 

■ 9m-2 

Go = 

9m—2 

9m —1 

90 ■ 

9 m—3 


9i 

92 

93 ■ 

90 
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An (m X m) matrix of the type Go is called a circulant matrix of order 
m or simply a circulant matrix. 

It is well-known that the generator matrices of QC codes can be trans¬ 
formed into blocks of circulant matrices by a suitable permutation of columns. 
Therefore, any 1-generator QC code is permutation equivalent to a linear 
code generated by a matrix of the form 

\ Gi G 2 ... Gi] 

where each G^ is a circulant matrix of the form (3). 

2.3 A New QC Z 4 Code 

Based on the search results for new codes over fields and over Z 4 , it is 
natural to search for new quaternary codes in the class of QC codes over 
Z 4 . Although this kind of search is carried out in MM and |23], a more 
complete search is still possible. For each odd integer m up to length 63, 
we produced all cyclic codes i.e., their generators, p(x),(free or non-free), 
based on the results described in section II-A. We then searched for new 
QC codes of the form {p{x),p{x) fi{x),... ,p{x)fi-i). In most cases we used 
/ = 2 (in a few cases we also let / = 3,4). Our search revealed a good (and 
new) quaternary QC code with parameters [ 86 ,4^® 2 ^, 55], whose Gray image 
(which is non-linear) is a binary (172,2^^, 55)-code. The best known binary 
linear code of length 172, and dimension 30 has minimum distance 54. The 
generators and the Lee weight enumerator of this code are as follows: 
g{x) = x^^ + 3x^^ -|- 2x^^ + 3x^^ -|- 2x® -|- 2x® -|- 2x^ + 2x® -|- -|- 2x^ -|- x -k 3, 

/(x) = x^®-|-x^’^-|-3x^®-|-2x^®-|-x^^-|-2x^^-|-3x^^-|-x^°-|-3x^®-|-2x^®-|-x^’^-|-x^®-t- 
2x^®-|-3x^'^-|-2x^^-kx^^-kx^^-|-2x^°-|-3x®-|-x®-k 3x^-1-2x®-|-x^-|-2x^-|-3x^-kx-k 1 
h{x) = 1 so that g{x)f{x)h{x) = x^^ — 1. 

Let p{x) = f{x)h{x) -k 2/(x), then the polynomial p{x) generates a free 
quaternary cyclic code with parameters [43,4^®2^, 16]. The search revealed 
that with the choice of fi = 2x^^ -k x^^ -k x^° -k 2x® -k 3x^ -k x'^ -k 3x® -k 3x® -k 
3x^ -k 2x^ -k X, the 1-generator QC code generated by (p(x),p(x)/i(x)) has 
parameters [86,4^®2^, 55]. Its Lee weight enumerator is given below where 
the bases are weights, and exponents are number of codewords of that weight 

0155774^015915'^36985g52895gl3244 00246390^436020274691 
03132870 0^233877 053741OO006141690'^97O8540gl5O2291 
002252598 y0332O2O2 '^^4791318 y26689811 739186262 74 I 2274866 
7515998236 7020463442 77255984 I 6 7g31106974 7g36948696 g043O8O625 
3^48872424 g254121520g358775152 g4622578513064430426 g065299285 
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§-^64550138 gg62322437 gg58728454 gg54154888 48850752 g242923718 

g337050520 g431176720 g525516630 gg20478707 g';7l6029368 ggl2290346 
gg9187466 ]^gQ6707312 ]^g]^4753392 ^023279137 ^032255178 ]^04498636 
Iq^982292 200634379 ^ 07382872 2 03227341 20gl3459O 2 2076067 
2224145222221930213105782246665225344022511182171032223172 
120129121 ® 6122®®1292 

3 Quaternary Codes from Inverse Gray Map 

The Gray map is usually used to obtain binary codes (usually non-linear) 
from quaternary codes (usually linear). However, we can also use its inverse 
to obtain quaternary codes (most likely non-linear) from a given binary code. 
If we take a binary code with parameters [2n, 2k, d] then the inverse Gray 
map yields a quaternary code (which is most likely to be non-linear) with 
parameters (n,4^,d). Taking advantage of existing databases for binary 
linear codes, we considered quaternary codes obtained this way from best 
known binary codes. This method contributed thousands of (non-linear) 
codes to the database. 

4 Genetic Code Search 

It is well known that computing minimum distance of an arbitrary linear 
code is an NP-hard problem |2S]- This result gives an insight about why 
there does not exist an efficient, general purpose search algorithm to find 
good linear codes. All known search methods/algorithms for linear codes 
work well in some special cases. Recently, a new search algorithm has been 
introduced that works well for most parameter ranges over small fields [4] . In 
a large number of cases the algorithm produced linear codes with best known 
parameters, and in several cases generated new codes (“record breakers”). 
In our work we adopted this algorithm for the ring Z 4 and the Lee metric. 
We refer the reader to [4] for further details. Originally, as implemented for 
the field case, PDG did not work very well for Z 4 . Therefore, we introduced 
changes inspired by genetic algorithms. We start off with an empty 
code and gradnally expand it. At each step we examine multiple 
mutations of a single generator matrix. Thus GCS is not a typical 
genetic algorithm, in the sense that it operates on a single element 
and crossover between generator matrices is not considered. Here 
we present the details of GGS for free codes, i.e. K = Ki,K 2 = 0. 
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Initialize Set the input parameters and initialize 
sets and variables. 

Set N,K,T; 

Use the binary record table to determine D] 
BitShifts = {1, 2, 3}; 

^ = {}; 

G = [0];t = l;fc = 1;_ 


The length N and dimension K are the two input parameters. Based 
on those we can determine a desired Lee minimum distance D 
using the record table for binary codes. In addition, we also specify 
a small integer T, 1 < T < N — K. Greater T implies a better chance for 
the construction of a good code but that benefit comes at the expense of 
increased computational time. The output of the algorithm is a linear 
code over Z4 of length N and minimum Lee distance d > D. 


General GCS Algorith 

Initialize 

while {{t < T) and {k < K)) do 
S = S\J{{K + kN), ...,{{k + l)N- 1)}; 
Gold = [ 0 ]; Gnew = [ 0 ]; 

Gtemp = [0]; GoZd [/?][/?] = 1; 

d= 1; 

Search for a suitable matrix Gold 
if {d > D) then 
G = G + Gold] 
k = k + 1; 

end if 
end while 


The search for a suitable mutation matrix Gold presents the heaviest compu¬ 
tational task. As the size of the set S increases with successive dimensions, 
so does the number of possible mutation matrices. This is a key differ¬ 
ence between PDG and GCS. The process terminates either when a suitable 
mutation matrix is found or the specified level of T is reached. 



Search for a suitable matrix Gold 
while ((t < T) and {d < D)) do 
increment Jt = true; 

Redundancy Shifts = The set of all subsequences 
of BitShifts of length t; 

Positions = The set of all subsets of S of size t; 
for p in Positions do 
Gnew — Goldj 

for r in Redundancy Shifts do 

for i = 1 to t do 

Gnew\p{i) ^ N][{p{i) mod iV) + 1 ] = 

= Goidlpii) N][{p{i) mod iV) + 1] + r{i); 

end for 

Gi = G + Gnew] G =< Gi >; 
if {MinimumLeeWeight(G) > d) then 
d = MinimumLeeWeight{C); 

Gtemp — Gnewi 

increments = false; 

break p; 
end if 
end for 
end for 

Gold — Gtempj 

if increments then 
t = t 1; 

end if 
end while 

Below is a table of small quaternary codes obtained with GCS whose Lee 
distances are equal to the minimum Hamming distances of the corresponding 
binary linear codes. We call such a code “decent”. This is significant, 
considering that all binary linear codes up to length 32 are optimal. Based 
on our experience and the results from the literature, constructing decent 
codes is a challenging task. Besides, many of the decent Z 4 codes may very 
well be regarded new. 
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Free decent quaternary linear codes. 

N 

Ki 


10 : 

1,2,3,4,5, 

6 , 7, 8 , 9 

11 : 

1,3,4,7,8, 

9,10 

12 : 

1,2,3,4,5, 

8,9,10,11 

13 : 

1,2,3,4,5, 

7,8,9,10,11,12, 

14 : 

1,2,3, 6 , 8 , 

10,11,12,13 

15 : 

1,2,3,4, 7, 

9,11,12,13,14 

16 1 

,2, 6 , 7, 8 , 9 

17 : 

1,2,3,4,5 


18 : 

1,2,3 


19 1 

,2,3 


20 : 

1,2,3,4 


21 : 

1,2,3,11 


27 : 

1,2,3 



5 Conclusion 

In this work, we introduce a new database of Z 4 codes that is available 
online that can be conveniently updated by researchers. The database has 
been populated using several different search methods. We present a survey 
of some of the recent and promising methods to find new quaternary codes. 
Search with one of these methods has yielded a good quaternary code. We 
also introduce a new search method that has yielded decent codes. We 
invite researches to search for new quaternary codes using known methods or 
devising new ones, and update the database with any new codes discovered. 
There is much room for improvement on this database. 
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